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Abstract. We briefly review the most relevant aspects of complete integra- 
bility for classical systems and identify those aspects vfhich should be present 
in a definition of quantum integrability. 

We show that a naive extension of classical concepts to the quantum frame- 
work would not work because all infinite dimensional Hilbert spaces are uni- 
tarily isomorphic and, as a consequence, it would not be easy to define degrees 
of freedom. We argue that a geometrical formulation of quantum mechanics 
might provide a way out. 

PACS:03.65.Ca, 02.40.Yy, 02.30.Ik 

1. Introduction 

In classical Mechanics, completely integrable systems have been widely studied 
in the past thirty years, many aspects have been considered and several approaches 
have been proposed ( IHH [H] ) . Many soliton equations have been shown to 
admit a Hamiltonian formulation and to be completely integrable. Many completely 
integrable systems have been shown to arise as reduction of "simple" systems, very 
often associated with geodetical motions on Lie groups ([IS]). In this respect, they 
have been shown to be associated with Lie-Scheflers systems [35J[n], i.e. systems 
whose solutions admit a superposition rule [5]. Thus we may safely say that clas- 
sical completely integrable systems are well understood. On the other hand, it is 
commonly accepted that our description of the physical world should rely on Quan- 
tum Mechanics. As Dirac puts it: "classical mechanics must be a suitable limit of 
quantum mechanics" . 

Inspired by this statement, it is quite natural to look for a definition of quantum 
integrability whose appropriate limit would reproduce the classical definition and 
the multifacets situations one knows for classical systems, namely: multiHamil- 
tonian descriptions, recursion operators. Lax pairs and so on [24l[34l[25]. On the 
mathematical side, the quantization of soliton models has given rise to the beautiful 
structure of Yang-Baxter relations and the related formulation of quantum groups. 

However, there are several problems if one tries to implement in a naive way the 
classical situation. Roughly speaking, the problems can be reduced to the fact that 
in infinite dimension all Hilbert spaces are isometrically isomorphic and therefore 
all C*-algebras of bounded operators acting on them are isomorphic. In particular, 
the idea of using the maximal set of commuting operators to imitate a maximal set 
of commuting first integrals of the classical situation does not work. The reason 
is a theorem by von Neumann stating that any two commuting operators can be 
written as functions of a third one, i.e. out of the finite dimensional maximal set of 
commuting operators we can find just one operator allowing us to express all the 
others as functions of it f[55]l. 
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From this point of view, what is missing is a way to distinguish, say, the Hilbert 
space of states of one particle from the Hilbert space of the states of many particles. 
In other terms: in the quantum setting it is not clear how to define the notion of 
degrees of freedom. Kirillov [29 , for instance, suggested the idea of introducing 
the notion of "functional dimension" . In some sense, we need to implement at the 
quantum level the notion of "functionally independent first integrals" . 

Besides these problems, other problems arise from the need of a sound definition 
of "suitable classical limit", i.e. to attach a mathematically reasonable meaning to 
the physicist's limit H ^ |43j . 

At the moment, the way we consider this limit depends on the quantum picture 
we deal with. For instance, the Schrodinger picture versus the Heisenberg one, 
the first dealing with equations of motion on a Hilbert space the second dealing 
with equations of motion on a C*-algebra. The two pictures are usually connected 
by the so called GiVS'-construction [17 . At the classical level the first picture 
corresponds to describe the dynamics as a symplectic vector field on a symplectic 
manifold, while the second one would correspond to the description of the dynamics 
as a Poisson derivation of a Poisson algebra. 

By using this analogy we may consider the Hermitian structure on the vec- 
tor space of states as the mathematical structure corresponding to the classical 
symplectic structure and similarly the Lie- Jordan algebra structure on observables 
would replace the classical Poisson structure. 

However, before entering more deeply into the analysis of these analogies, let us 
consider what would be a reasonable approach to the integrability of a quantum 
evolution equation. Given the operator H in some concrete realization, say the 
Schrodinger picture, we look for a Lie algebra q and a group G that are "naturally" 
associated with H. Typically, this would mean that there is a unitary representation 
of G on the Hilbert space (the carrier space of our system) such that H can be 
identified with the action of an element in the universal enveloping algebra U{q) 
of the Lie algebra g. The structure of the unitary representations of G and their 
decomposition into irreducible ones then yields information about the spectrum 
and the eigenfunctions of H. This approach relies on the implicit assumption that 
the theory of unitary representations of G is sufficiently developed so that it may be 
used as a non-commutative Fourier analysis (non-commutative Harmonic analysis 
on G). Nowadays a lot is known for a large class of connected finite-dimensional 
Lie groups 

Fundamental tools for this analysis are provided by Mackey's imprimitivity the- 
ory, a generalization of techniques elaborated by Weyl and von Neumann within 
quantum mechanics characterizing such representations; and the Kirillov-Konstant 
correspondence between representations and coadjoint orbits. For the use of Lie- 
algebraic methods in the study of representations it is instrumental also the use 
of the space of smooth vectors (this originated with coherent states introduced by 
Schrodinger and further elaborated by Bargmann and Nelson); which is a module 
for the associative algebra [/(g). 

The connection with biHamiltonian classical systems arises with the introduc- 
tion of biunitary operators. They are associated with one-parameter groups of 
transformations which are unitary for alternative Hermitian structures on the same 
topological vector space of states (UHl HI]). The geometrical formulation of Quan- 
tum Mechanics turns out to be very useful to identify the relevant mathematical 
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structures so that we may, more easily formulate what we mean by alternative 
structures compatible with a quantum evolution. 

Let us recall briefly the fundamental aspects of classical integrability. 

2. Composition of model dynamics 

2.1. Classical integrability. From a classical point of view, the theory of inte- 
grable systems is very well established. In this section we shall just discuss a few 
of the best known properties which will become later important when compared to 
the analogue situation in the quantum framework. 

Let us begin considering the usual definition of a classical integrable system. 
Consider a classical system with n degrees of freedom, defined on some differential 
manifold M. For simplicity, let us consider a dynamics F defined via a Hamiltonian 
function H on the cotangent bundle T*M. We will consider the canonical Pois- 
son structure associated to the canonical symplectic form, and the corresponding 
Poisson bracket defined on C^{T*M). 

Definition 1. The system is said to be integrable if there exists 

• a set of n functions where fk G C°°{T*M), which are constants 
of the motion 

= {H,fk} Vfc 

• functionally independent dfi A • • • A c?/„ ^ 

• which are in involution, i.e. 

{/fe,/,} = Vj,fc 

Of course we could consider more general situations to describe systems whose 
dynamics can be explicitly integrated. In principle, just the fiow of the system 
on the given manifold is to be considered as the object of interest. The geomet- 
rical structures around it may be changed and actually this leads to important 
consequences in some situations, as in the case of biHamiltonian systems. 

There a few properties of these systems which we would like to consider more 
closely. First of all, one of the most relevant properties is that the dynamics for 
such a system is always integrable by quadratures. For instance, we can write the 
dynamics in a trivial way in a new set of coordinates, defined by the set of the 
constants of the motion and their corresponding conjugate variables. These can 
always be found, at least locally [3| [T| |4T|: 

Theorem 1. There exists a set of coordinates {Ij,4'j}^=i such that the mapping 
: {q^,Pi) I— > {IhTfj^k) is a local symplectomorphism and such that the vector field 
$,(r) is written as 

^*ir)^J2^k{i)^ (1) 

k 

This is one of the relevant features we would like to be able to recover at the 
quantum level: the possibility of integrating completely the dynamics in a construc- 
tive way. This is one of the more involved points in what regards the definition of 
quantum integrability, because, within the quantum setting, it requires nonlinear 
transformations of non commuting variables, and we shall discuss it briefly. 

There is another aspect of classical integrability which we would like also to 
mention here. The definition above is somehow restrictive from certain points 
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of view. If what we want is just to integrate the dynamics, we could drop the 
requirement that $ be a symplectomorphism and just require that ([l]) is a vector 
field which is nilpotent of index two. This would allow to cover much more general 
situations. Let us review this quickly: 

2.2. The meaning of Arnold-Liouville theorem. From our point of view, 
Arnold-Liouville's theorem (AL in the following) can also be read in a different 
way. One should interpret it as a way of relating a given system, defined in the 
original manifold M with a set of systems defined on suitable associated model 
manifolds {-/Vj} via "projections" <i>j : M — > Nj, defined in action-angle variables, 
whose dynamics is trivially integrable. The role of the n constants of the motion 
in involution is just to give rise to the transformation which allow to define these 
"projections" . This transformation is nothing but the set of mappings defining each 
torus as a leaf of the foliation of M according to the level sets of the momentum 
mapping. Each leaf may be therefore identified with our manifold Nj, and the 
transformation 



is the corresponding mapping <i>j. 

The aspect that we want to point out is that, from this perspective, AL construc- 
tion is analogous to the transformation of a linear system into Jordan blocks. The 
submanifolds characterized by $ and v are now the analogues of the Jordan blocks. 
There is an important difference however: AL theorem shows that by using non- 
linear canonical transformations we can relate generically any complete integrable 
systems with two types of models: 

• "Harmonic oscillator" -type for those systems whose dynamics is completely 
bounded 

• "free particles" -type for those systems with unbounded motions. 

The decomposition of the dynamical vector field is also obtained as above. Sup- 
pose that we are working with action-angle variables on the classical manifold M 
that we assume to be symplectic. Classical dynamics is described in terms of a 
vector field, which is the derivation associated to the classical Hamiltonian H. 
Suppose we have an integrable system, whose Hamiltonian vector field is written, 
in the usual way, in action-angle variables {/j, 0j}j=i...„ as: 



where each Vj e M is the corresponding "frequency" and we assume ipj € (— 7r,7r) 
(we assume hence each single dynamics will be bounded) . This vector field describes 
the motion of a "classical Harmonic oscillator" and its integral curves are contained 
in the corresponding toroidal manifold which serves as level set for the momentum 
mapping. On any torus, the dynamical vector field has the same form, since the 
frequencies are constant, however, in general, these constants may depend on the 
torus 



We should point out that by means of a nonlinear transformation we are going 
from a system described, say, by the matrix 



{fe,p')}^(/„0,) 



Remark 
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in the {q,p) coordinates to a system described by the nilpotent matrix 

1 


in the (<I>, z/) coordinates. Thus wc arc considering a nonhnear transformation which 
relates two linear dynamical systems described by completely different matrices. 

We may also describe the structure above in a generalized fashion by requiring 
the existence of a set of functionally independent functions {Ji,--- ,J2n} which 
satisfy 

CriCrJj) = 0, 

with r being the dynamical vector field. Wc obtain in this way a description in 
terms of generalized "action-angle functions", which need not be pairwise "con- 
jugate", instead of the local "action-angle coordinates". It means, in particular, 
that a completely integrable system is nilpotent of order two, on some set of func- 
tionally independent functions which separate points of the carrier space. In these 
"coordinate functions" the formal exponentiation gives 

e^^^Jj = Jj+tLrJj, 

i.e. the motion is linear in time. If we replace Jj with we recover the 

dynamics on each torus: 

Let us consider more closely integrable systems which are not described by 
isoperiodic motions. The vector field on each torus is always of the same type 
as above, but now the frequencies exhibit a dependence on the action variables, i.e. 
j/j = • • • , /„). On each torus the vector field is "equivalent" to the vector field 

of the Harmonic oscillator, i.e. a linear combination of the angular vector fields, 
but the coefficients are different for each torus. As a consequence, wc can say that 
a classical integrable system is made up by a suitable combination of vector fields 
corresponding to Harmonic oscillators, where suitable means that the oscillator 
vector fields from each toroidal leaf arc "composed" by using "coupling constant" 
which are functions of the action variables. 

2.3. Conjugate systems and classical integrability. In this section we want 
to elaborate further, with some additional details and from a more geometrical ap- 
proach, the interpretation of AL theorem that wc just described. The idea is that 
the integrability of a vector field X (the dynamical vector field of our system) may 
be related with that of a set of models {Yj}, related to a transversal decomposition 
of the original manifold, whose evolution is well known. So we say that the inte- 
grability of each one of the Yj ensures the integrability of the vector field X, which 
is supposed to be decomposed in an analogous manner to the canonical "Jordan 
form" of a linear system. We will see in the following that this procedure offers 
many advantages. 

A crucial part of the present work concerns thus the relation between structures 

belonging to different spaces, particularly vector fields and algebras of functions. Wc 
will begin with the differential geometric preliminaries needed for the understanding 
of our considerations. 

Let us consider two differential manifolds, M and TV and a differentiablc mapping 
between them ^ : M ^ N We suppose that we have a vector field on each manifold 
(the dynamical vector field for instance), and we want that $ relates the different 
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set of orbits on each space, in a 1:1 correspondence, without requiring that they 
have the same parametrization. The general expression for this correspondence is 
summarized in the expression: 

T$ o X ^ /r o $ (2) 

where T<i> is the induced tangent bundle map acting on vector fields, i.e. sections 
of TM and TN respectively, and / is an arbitrary function on N. In other words, X 
and Y are (<&, /)-related, i.e. X and fY are <i>-related (this is a way to formalize the 
statement that the orbits are in one-to-one correspondence up to parametrization) . 

When this mapping is such that / = 1 we have a 1:1 correspondence not only 
for orbits, but also for integral curves, i.e. the image by 3> of any integral curve of 
X on AI is an integral curve of Y on N; in this case we say simply that X and Y 
are <i>-related. 

We will use now this kind of mappings, suitably combined, to decompose the 
dynamics on a given manifold, defined through a Hamiltonian function or a Hamil- 
tonian vector field, in a set of well known dynamics on special manifolds which we 
will call "models" . 

In fact we will not use the dynamics on the full manifold but just on an open 
dense submanifold, extending suitably. 

Let us apply now some general considerations. Let M be a symplectic manifold 
(it is not necessary, but it will be so for our applications) and let X E X{M) be the 
dynamical vector field. Let us call model systems A4j, some pairs {Yj,Nj) where 
Nj is a differential manifold and Yj G X{Nj). We suppose we can define a covering 
{Mj} of AI along with a set of diffeomorphisms: 

$j : A4j CM ^ Nj 

which posses property 

Now we want to use these sets Nj and these vector fields Yj to describe the 
dynamics corresponding to X. We can choose two different ways for doing it: 
working at the level of vector fields and manifolds or working at the level of functions 
and their algebraic structures: 

a) We can work at the level of vector fields, where we know the dynamics is 
described by X, and we can try to find a suitable set of vector fields Yj in 
such a way that we can integrate this dynamics with their aid. We suppose 
that (j2]) holds for any Nj , and that a decomposition may be chosen in such 
a way that we could write: 

(3) 

i 

where Cy^hj — Q "^i ^ j and Yj e X(M) is the vector field on M <I>j-related 
with Yj. 

In this set, we define: 

Definition 2. We say that the model set M. = {{Nj, Yj)} gives a transver- 
sal decomposition for X(M) when we can write: 



X(M)(m) = 0[x](iV,)] (m) 



TOWARDS A DEFINITION OF QUANTUM INTEGRABILITY 



7 



which means that at any point of the manifold M the "lifts" of the vectors 
fields from Nj to M define a basis of the corresponding tangent space T^M. 

Thinking in terms of the C°° (M)-m{Klulc structure of the set of vector 
fields of M, it is also interesting to further analyze this decomposition. We 
could say that the set of vector fields {Yj}, the lifted vector fields, defines 
a system of generators for the module X(M) with coefficients in C°°(M). 

When we work in classical dynamics, the physical observablcs are func- 
tions C°°(M), and therefore we want to describe the previous construction 
in terms of functions. We can restate the previous requirements in the form 

£j.$*(/) = $*(Ly^./) V/eC°°(iV,). 

b) We can also work at the level of functions on each manifold with the algebra 
structure provided by the Poisson bracket. The dynamics now is described 
with the Hamiltonian function and its adjoint action inside the Poisson 
bracket: 

The set of models now is just specified by {Nj} (and the corresponding 
mappings since we write the definition of transversal decomposition 
as: 

Definition 3. Given M and X as above, and the set of models and map- 
pings {Nj,^j) we say that they define a transversal decomposition of an 
open dense submanifold of M when the set ®jd^*(C°^'(Nj)) generates 
X*{M) as a module (the set of 1-forms on M). As above we are requiring 
that 

d^*{C^{Nj)) n rf$*(C°°(7Vi)) = , i ^ i . 

And moreover, we say that a vector field X g 3£(M) is decomposable 
according to the set of models if 

i.e. each subring <^*{C°°{Nj)) of functions is stable under the action of the 
dynamical vector field. 

It is important to remark that both approaches are not trivially related, 
since the set of Hamiltonian vector fields (which is the set arising in this 
latter one) does not define a C°°(M)-module, and therefore there is no 
possible transversal decomposition associated with them. 
Remark 

From this point of view, even non-Hamiltonian systems may be considered in- 
tegrable. In the first cases we have discussed objects corresponding to the infinite 
dimensional Lie algebra of Hamiltonian vector fields. But it is also possible to 
formulate things directly at the level of the C°°(M)-module of vector fields as a 
whole. Suppose that the dynamical field can be written as a linear combination of 
derivations that are inner when restricted to a torus, with coefficients depending on 
the action variables. Using the preceding formalism we say that in the combination 
above we do not need to impose that the functions hj belong to ^j{Nj) but just 
that: 

Cyh = (4) 
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and therefore the derivations corresponding to X need not be inner (corresponding 
to a Hamiltonian on M) on the whole manifold, but that it coincides on each torus 
with an inner derivation which may depend on the torus. We could have a function 
hj = hj{Ii, ■ ■ ■ , /„) which is constant on each torus, but that depends on all of the 
constants of the motion, the equation Q above implies only that the Hamiltonian 
commutes with the function hj and therefore its most general form is an arbitrary 
function of the constants of the motion of the system. 

To put it differently, we could use one-forms which are not closed to describe 
vector fields by means of the symplectic structure. What we require is only that 
the one-form is expressed in terms of action-angle variables. In different terms, in 
the equations 

we may drop the requirement that are partial derivatives of a function (the 
Hamiltonian function) [5]. 

2.4. Examples of classical integrable systems. Now we are able to give a few 
examples which illustrate the situation in the classical case. We can consider two 
families: 

(1) Systems whose orbits are on tori. Model systems in this case would be just 
combinations of one dimensional Harmonic oscillator. Examples of this 
group would be: 

• The n dimensional Harmonic oscillator 

• The q-deformed oscillator, which would be the simplest non trivial 
example. 

• The closed orbits of the Kepler problem. 

(2) Systems whose orbits are on cylinders (i.e. on K*^ x T'^-^y xhe model 
systems now have to combine the one dimensional Harmonic oscillator and 
the free particle, to include the unbounded motion. As examples of this 
case we could mention: 

• The unbounded orbits of the Kepler problem. 

• The Calogero-Moser type potentials. 

2.4.1. Example: the Harmonic oscillator. Let us develop now some of these exam- 
ples. Let us begin with systems defined on tori. 

For this kind of systems, the preceding description of Liouville integrability is 
straightforwardly adapted. Let us take as a first example the n dimensional Har- 
monic oscillator: 

k 

If we select as first integrals the partial energies, we can define: 

I J = hj = ^{p^j+Ujq]) 
and the Hamiltonian is therefore written: 

H = Y.h, = j2h 

k k 
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Liouvillc- Arnold theorem provides us also with the angle variables, and the con- 
formal factor F]^ is obtained following the expression above as: 

/.-^=-. (5) 

Obviously, the Harmonic oscillator turns out to be integrable, and its dynamics 
linear. 

The following example is the first non trivial modification: the q-deformed clas- 
sical oscillators that can be found for instance in [39 . The Hamiltonian may be 
written as: 

sinh hakal 



k 



sinh h 



where ak^^{^'^ + V^Qk) ■ 

The invariant now is the classical "number operator" , but the conformal factor 
will not be a numerical constant, but a constant of the motion, i.e. a function of 
the invariants: 

h = oikOi*k 

Jk = -WT ^ ^ ■ 1 fc coshtofcafe 
oik smh h 

The role played by this conformal factor is quite clear, it represents a frequency 
that is not constant on the phase space, but it is constant on each manifold /fc, 0^. 

2.4.2. Conformal factors. On M^, with coordinates {x,p) we consider the dynamical 
system: 

x = pf{x^ +p^) p^-xf{x^+p^) 

It is immediate that the orbits of this system are all circles with center at the 
origin. We can integrate our system for any initial condition (xo,po): 



x{t) 
Pit) 



exp 



tfi4+pl) 



1 




xa 
Po 



It is clear that we are using the Harmonic oscillator system 

X ~ Lup 

p = —LUX 

to find our flow and replace the constant frequency uj with a constant of the motion 
/. If, however, we look for a diffeomorphism "conjugating" our starting system 
with the Harmonic oscillator, we would find that no such a diffeomorphism can exist 
because we know that the period of each orbit is invariant under diffeomorphisms, 
and our two systems have different periods for each value of the radius. Therefore 
the conformal factor will help in declaring our system "conformally equivalent" to 
the Harmonic oscillator. Thus without the conformal factor, each / would provide 
a different equivalence class. 
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2.4.3. The Kepler problem. We consider now the case of the Kepler problem. In 
principle the idea is simple, we have just to proceed as we did above. As a first 
point, the action-angle variable have to be determined. Following [IJ we find that: 




^3 = P^ (8) 
, _ ^-m^k^r^ + 2mfc( Jl + J2 + Jgj^r - ( Ji + J2 + J3)2( J2 + JsY 



( Jl + J2 + J3) 



2 



mkr — ( Jl + J2 + J 



3^ 



^2 



(Jl + J2 + J3)ViJl+.n + J3Y -iJ2 + J3Y 



(9) 



. [TOfcr-(J2 + J3)2](Ji + J2 + J3) . (^2 + ^3) cos ^ 

arcsm = arcsm — (10 ) 

J^oote , ^ ^^^^ 



ViJ2 + JsY ~ Ji 

Analogously, we can write the Hamiltonian in terms of these variables: 

""^"(XtItj^ ^''^ 

From our previous discussion, we can compute now the conformal factor for this 
case, simply as the derivative: 

" dJk 

For this case there is just one conformal factor: 

(Jl + J2 + J3)'^ 
and the related vector field acquires the form: 

2mP / d d d \ 

^ {Ji + J2 + JiY\d<t^i^ 04^2^ dh) ^ ' 

Now we have to consider how to extend these structures to the quantum domain. 
It is clear that the geometrical formulation of Quantum Mechanics allowing for 
nonlinear transformations on the space of pure states, the space of rays or the 
complex projective space, will be more appropriate to carry on our considerations. 

3. The problems of Quantum integrability 

Let us try to adapt to the quantum setting the results of classical integrability 
we have analyzed. We recall that in order to be able to explicitly integrate the 
system, the main properties and requirements of AL theorem were 

• the number of degrees of freedom of the system, which determines the 
dimension of the manifold where the system evolves, is equal to the number 
of functionally independent constants of the motion in involution. 
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• it is important to characterize the geometrical objects the system may be 
associated with. The reason for this is that the existence of different ge- 
ometrical structures compatible with a given dynamics may provide new 
constants of the motion (as it happens with classical biHamiltonian sys- 
tems, for instance) 

• the notion of independence of the constants of the motion is crucial for them 
to be used as coordinates in the action-angle description (and therefore they 
must be functionally independent). 

All these properties, well known and controlled in the classical framework, are 
not without ambiguities when transferred to the quantum setting. Let us study 
them in some detail: 

3.1. The notion of degrees of freedom. Let us now turn our attention towards 
the von Neumann theorem we mentioned in the introduction. It reads as follows: 

Theorem 2. Given a Hilbert space Ti. and two observables A, B, if they commute, 
there exists a third operator C such that both A and B are Junctions of C. 

Thus if we look for observables which commute among themselves and with the 
Hamiltonian, as the naive generalization of AL theorem would suggest, we face 
the problem of dealing with their independence since the theorem we have quoted 
asserts that they will be functionally dependent on just one operator C. We can 
not then relay on the number of degrees of freedom of the system, for this is a 
concept which does not have a clear meaning in the usual quantum setting. Let us 
elaborate a bit in this direction. 

We are going to see how it is not possible to define a reasonable notion of degree 
of freedom in the case of quantum systems. To this aim, we shall consider an 
example taken from a paper by Weigert [56 and further elaborated by Hietarinta 

m- 

Consider the case of a one-dimensional Harmonic oscillator 

H ^ huj ^fltfl = huj + . (15) 

We recall that there is a countable set of eigenstates for the number operator h: 

h\n)=n\n) n G No {0, 1, 2, • • • } 

which constitute an orthonormal basis for the carrier Hilbert space of states of our 
system. For the creation and annihilation operators we have, in this basis, 

a\n) — \/n\n — 1) a^\n) = + l\n + 1), 

along with 




We would expect that this system, being one dimensional, should have only one 
constant of the motion. However, we are going to show how it is possible to define 
two constant of the motion in involution by using the following trick. 

Let us replace the label n by two indices ni and 712, both taking all non-negative 
integer values and related to n by the following relation 

n = ni + i(7ii + n2)(ni + 712 + 1). 
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We can list then a few correspondences: 

n = — > (ni,n2) = (0,0) n = l — ^ (ui, Us) = (0, 1) n = 2 — > (m, ria) = (1, 0) 

n = 3 ^ (m, na) = (0, 2) n = 4 ^ (ni,n2) = (1, 1) 

Considering a two dimensional grid (in an analogous way we can repeat the 
argument with three operators getting thus a cubic grid, etc) we find the following 
diagram 



Every pair (ni,n2) determines exactly one integer number n and vice- versa. 
Therefore it is possible to label the states \n) by using the two integer- valued labels 
(ni,n2): 

|n) = |ni,n2). 
We can find, obviously, an orthonormal basis 

along with a decomposition of the identity 

|ni,n2)(ni,n2| = ^ = I 



These states can be considered as the eigenstates of two operators, fii and n2, 

independent and hence commuting among themselves. 

The eigenvalues of the Hamiltonian operator may be written as 



En = hu) ( ni + ^{ui + n2)(ni + ^2 + 1) + ^ 



ni,n2 e No 



The states |0), |1), |3), • • • ||fc(fc + 1)) belong to rii = and therefore the eigen- 
value has infinite degeneracy. Similarly for fi2, for which we can write: 



n2 = Y.Y.^P (l+\{k + l){k + l + l\ . 

1=0 k ^ ^ 



Here P{m) = \m){m\, i.e. this represents the operator giving the spectral decom- 
position for h2- 
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It is clear that [ni , •0,2] = and that the total Hamiltonian operator can be 
written as 



H = hu{hi + ^(ni + n2)(ni + na + 1) + ^ ) . (16) 



Obviously, 



[H,ni]=0 = [H,h2]. 



But this means that 



Lemma 1. The one dimensional Harmonic oscillator above has two constants of 
the motion commuting with each other. 

Wc can consider the eigenstatcs \n) as belonging to the Hilbcrt space which is 
the tensor product Tii(^Ti.2, Ti-i being the Fock space corresponding to each number 
operator. Both operators h\ and n2 have infinite degeneracy, but the states of the 
composed operator h has none. 

H\k) = k\k) fee No = |ni) ® [r?,2> ni,n2 € A^o- 

This has been done for the particular case of the Harmonic oscillator, but as we 
are going to sec, any one-dimensional quantum system (with discrete spectrum and 
eigenvalues bounded from below) can be given previous formulation. Consider for 
instance a system described by the Hamiltonian 

2s; + ''(')^ 

and assume its spectrum is lower bounded as i?o < ^'i < -£'2 < • • • . Denote the 
elements of the orthonormal basis of eigenstates as {|i?j)}j=o,i, -. 

We can define now linear operators on the complex Hilbert space spanned by 
the above vectors by setting: 

A\Ek) = Vk\Ek-{); A^\Ek) = Vk + l\Ek+i), 

along with 

Vn! 

It is now obvious that the product A = N has the same eigenvectors as the 
original Hamiltonian H and thus H can be written as a function of N . Therefore 
the same considerations we gave above for the Haxmonic oscillator can be given for 
a generic Hamiltonian (with the properties on the spectrum we discussed above). 
But this implies that the system admits two constants of the motion Ni and N2, 
in terms of which the Hamiltonian can also be written. 

More generally, given a maximal set of commuting operators, uniquely identifying 
an orthonormal basis for the Hilbert space, once we order these states by means of 
No, we can define an operator 

(") 

Thus, any element of the family of commuting operators will be an appropriate 

function of K. 

It is clear that these peculiarities are to be traced in the fact that all infinite 
dimensional Hilbert spaces are isomorphic, indeed unitarily isomorphic. If \n)(a) 
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(with n G N) provides an orthonormal basis for 7i(a) (where a ~ 1,2, • • • ), we can 
estabhsh this correspondence by setting \n)2 = T\n)i. That being so, any operator 
on 7Y(i) has a corresponding operator A(^2) = on 7i(2)- Thus also the 

corresponding C*-algebras are isomorphic. This means that, at an abstract level, 
we can not implement a notion of degree of freedom, because the Hilbert space of 
one particle is isomorphic with the Hilbert space of many particles (Balachandran: 
bringing up a quantum baby). 

With these examples we hope to have convinced the reader that the notion of 
degree of freedom in the framework of Quantum Mechanics does not have, at the 
moment, a sound mathematical definition. 

3.2. The integration of the dynamics. From the considerations we have made 
in Section |3.1| it is clear that we need a notion of functional independence for 
differentiable functions to be able to distinguish the carrier space of one particle 
from the carrier space of many particles. 

To the extent that the differential calculus is developed nowadays, to deal with 
(local) derivations we need commutative algebras. If we start with algebras of 
operators on some Hilbert space, the spectral theorem asserts that every Hermitian 
operator is unitarily equivalent to a multiplication function (operator). 

Thus, out of a family of commuting Hermitian operators we get a family of multi- 
plication functions. According to the Gelfand-Naimark theory, the point spectrum 
of these commutative functions, with some additional requirements, can be given 
the structure of a smooth manifold. This seems to be the main route to extract a 
"classical manifold" from a quantum carrier space. (As a matter of fact, the theory 
provides us with a Hausdorff topological space; to get a smooth manifold one has 
to select a subset of states, often C°° or analytic vectors) p!5l[5TlH5] . 

In some other approaches these aspects may be hidden in other assumptions and 
they appear in disguise. In many instances they appear very openly because the 
Hilbert spaces are required to be concretely realized as square integrable functions 
on some measure space, the algebra of operators is identified from the beginning 
with differential operators (in this case the commutative algebra of multiplications 
functions i.e. differential operators of degree zero, is immediately identified). A 
further choice of a specific quantum system, i.e. a specific Hamiltonian operator 
with a corresponding domain of self-adjointness, will select further topological and 
differential structures of the "classical manifold". 

In what follows we make few considerations on a procedure which subsumes most 
of the existing approaches to quantum integrability. We consider the equations of 
motion on some Hilbert space Ti., in the Schrodinger picture 

i^^V' = Hip, 

or with the equations of motion on the algebra of operators (Heisenberg picture): 

at 

We look for a Lie algebra g and a related Lie group G that are "naturally" 
associated with H. 

Typically, we look for a unitary representation of G in the Hilbert space Ti 
and require that the Hamiltonian operator can be identified with the action of an 
element in the universal enveloping algebra U{g), g being the Lie algebra of G. 



TOWARDS A DEFINITION OF QUANTUM INTEGRABILITY 



15 



The structure of the unitary representations of G and their decomposition in 
irreducible ones yields information about the spectrum and the eigenfunctions of 
H. Let us try to be more specific by recalling a few facts about representation 
theory: 

A unitary representation of G in the Hilbert space 7i is a homomorphism tt of 
G into the group of unitary transformations U (7i) with the following continuity 
condition: 

The map G x Ti ^ Ti defined as {g,'>p) TT(g)'ip) is continuous. The homomor- 
phism condition is equivalent to requiring that 7r(g) is a continuous linear operator 
on H for all g E G, that 7r(e) = I (e being the identity element of the group) and 
that the product is preserved, i.e. Tr{gig2) — 7r(gi)7r(g2) for all .gi,.g2 G G. 

If tt' is another representation of G in 7Y', a continuous linear map T : 7i — > 7i' 
is called an intertwining operator if T o 7r(g) = TT'{g) o T for all g G G. This 
would be the "quantum" version of our "related" vector fields considered in the 
classical setting. If T is an isomorphism, tt and tt' are said to be isomorphic. Given 
a representation tt, the continuous operators T : 7i — s- 7i which intertwine tt with 
itself, constitute an algebra, called the commuting algebra of the representation tt. 
For unitary representations, the isomorphic representations are unitarily equivalent. 

For a unitary representation tt : G — + U{Ti), we may consider the set of G°" 
vectors (denoted as 71°°) and the set of analytic vectors (denoted as TC^). According 
to Nelson [48], both of them are always dense subspaces of the Hilbert space TL. 
Thus it is possible to derive a representation of the Lie algebra g of G on these 
subspaces by setting 

TT : fl ^ End(7^°°); 

given by 

■h[x)il) = lim - (7r(exp(te))V' — "0) ■ 

This representation has a natural extension to a representation of the universal 
enveloping algebra U{q), which is still denoted by tt. 
If we define an involution on U (g) by setting 

X* = -X yX € 0, 

the extension tt to U{g) is a ^-representation in the sense that for all Y E U (g) 

^(y*)c(^(y))*. 

It is possible to introduce states on the universal enveloping algebra by means 
of linear functionals. A functional <i> : [/(g) ^ C is said to be a state if 

<i>{x*x) > yx e uio). 

This imphes that ^{X*) = ^{X) for all the elements X e U{g). Convex combi- 
nations of states are also states, and thus they form a cone which defines a partial 
ordering on the linear space of all self-adjoint linear functionals on U{g). Often, we 
consider normalized functionals defined as 4>(I) = 1. The states that can not be 
written as convex combinations are extremal states and are called pure states. 

By using the Gelfand-Naimark-Segal (GNS) construction, with any state we 
associate a Hilbert space H$, and a *-representation a of U{g) by (unbounded) 
operators on a dense subspace V C 7i$. There exists also a vector G I? such that 
for all X G [/(fl), 

{a{x)n\n) = <^>{x), 
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and V = (j{{U{Q)}).n. 

If $ is an analytic state, the representation is an irreducible unitary represen- 
tation of the algebra if and only if $ is a pure state. From the construction, it is 
clear that by acting with the group G on fi, we identify a submanifold of states in 
which may be pulled back to G and therefore put into correspondence with a 
Hilbert space contained in C^{G) (the Hilbert space of square integrable functions 
on the group with respect to the Haar measure). Thus, in conclusion, in this ap- 
proach a "classical manifold" appears again to be "built in" from the beginning 
via the smooth manifold represented by the Lie group G. Let us mention that the 
procedure considered in [21 shows very clearly that by using representations of G 
on spaces which carry alternative linear structures would give rise to "intertwin- 
ing maps" which are no more linear. In this way the "vector field" point of view 
would not be "optional" any more but would be mandatory. Moreover this would 
be a way to implement nonlinear transformations carrying from a Hilbert space 
description to another Hilbert space description not linearly related among them. 
A situation like this should be contemplated if we would like to have "alternative 
Poisson structures" non necessarily compatible among them. 

We shall mention now few examples where the previous constructions can be 
traced in the specific realizations of groups and algebras. 

3.3. Finite level systems and alternative descriptions. Within the geometri- 
cal formulation of Quantum Mechanics, a quantum system is described by a vector 
field r which preserves a Hermitian structure. The carrier space, when considered 
not as a complex vector space but as a real manifold, gives rise to a Riemannian 
structure, a symplectic structure and a connecting complex structure (these aspects 
will be further elaborated in Section 4). A vector field which describes quantum 
evolution has to preserve all of them: this means that it is not only symplectic, 
it is also a Killing vector field and generates a complex linear transformation. As 
the dynamics is linear and the structure above is translationally invariant, we can 
represent the previous compatibility conditions in terms of their representative ma- 
trices. 
We have 

Lemma 2. A matrix A is the matrix associated with some linear Hamiltonian 
vector field if and only if it can he decomposed into the product of a skew symmetric 
matrix A and a symmetric matrix H, i.e. 

A ^ A.H. 

Out of this decomposition we construct a Poisson bracket from A and the Hamil- 
tonian function from H. 

If we consider now any linear invertible transformation T : 7i — > 7i, we know 
the transformation law of the matrix A representing a linear transformation (A i—f 
TAT^^) and we can write, by considering A and H as representative of bilinear 
maps on account of the identification we are going to make 

At-. At^ TAT-^ = ThT\T^)-'^HT-^ = KtHt. (18) 

We have considered A to transform contravariantly and H to transform covariantly. 
If T is a symmetry transformation for A, i.e At — A, we conclude 

A = KtHt, 
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where, in general, Kt will be different from A. If A represents a Poisson structure, 
we would say that the transformation T is a symmetry for the dynamics which is 
not "canonical" . Thus the symmetry group for A, quotiented by the subgroup 
of canonical transformations for A would parametrize the family of alternative 
decompositions for A and consequently of the family of alternative Hamiltonian 
descriptions [47j . 

There are a few properties of these alternative descriptions which are immediate 
to list: 

• all odd powers of A are traceless, i.e. 

TrA^/c+i ^ 

and give rise to new Hamiltonian systems. 

• when A is generic, the commutant and the bicommutant coincide and are 
generated by all powers of A. 

• in the generic case 

e^-4'Ae-^(^")'=AA, 

for different values of A gives representative matrices of alternative Poisson 
structures. 

• in the quantum case, if T in addition commutes with the complex structure 
(i.e. it is the real form of a complex linear transformation), then it will 
also generate alternative Hermitian structures when H is positive definite. 
In particular, in the generic case, e^^ will generate alternative Hermitian 
structures (if we write A = iB where B is Hermitian, it is evident that the 
even powers of A are not skew-Hermitian) . 

By using collective coordinates, say {^j}, for the realification of the complex 
vector space, we have for the corresponding tensor fields: 

which explain the transformation properties we have used for A and H. We obtain 
thus alternative Hamiltonian descriptions. It is important to remark that the al- 
ternative Poisson structures are mutually compatible, i.e. any combination of them 
with real valued coefficients is still a Poisson structure. 

To be more definite, we shall consider a specific example given by a quantum 
system with a 2-dimensional Hilbert space (i.e. a two-level quantum system). The 
equations of the motion for such a system would be: 

S=[H,S] S= Siai + S2(J2 + Sacrg H = Ba-i 

The matrix may be considered to be the representative matrix of our dynam- 
ics, thus 




This case is not generic, and therefore we have to consider matrices in the com- 
mutant of which are skew Hermitian. Without going through the full procedure, 
it is not difficult to see that all alternative Hamiltonian descriptions may be written 
in the form 

h — a\ziZi + 02^2^2 fli, 0.2 G Hi; 
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and in real coordinates {zj = qj + ipj) we would have 

h = al{pl + qf) + al{pl + ql) 01,026 M; 

along with corresponding symplectic and Riemannian structures 

g — ai{dpi + dq\) + a\{dp\ + dgf) uj — afdpi A dqi + a\dp2 A dq2 

A further generalization arise if we let oi to become a function of ziZi and 02 a 
function of Z2Z2 (see [18]), or more generally 

h — a\{zizi, Z2Z2)ziZi + Z2Z2a\{ziZi, Z2Z2) ■ 

In this case the skew-symmetric part of the Hermitian tensor is still given by a 
symplectic structure, but the Riemannian part is no longer flat. If we use a Dar- 
boux chart for the new symplectic structure, we can define a new flat Riemannian 
tensor, and hence a new Hermitian structure, but it will not be connected to the 
original one by nonlinear "point transformations", i.e. transformations only on the 
"coordinates". To put it differently, the Darboux coordinates are adapted to the 
symplectic structure but are "not adapted" to the Riemannian tensor. The new 
Darboux chart allows to introduce a new linear structure such that the original dy- 
namics is still linear with respect to the new linear structure. On each symplectic 
plane the metric tensor is going to be 

{dp,f + pl{dcj,,f 

in the flat case, and there will be a conformal factor instead when aj is a function 
of the coordinates. 

By using the new linear structure associated with the Darboux chart, it is pos- 
sible to define a new realization of the unitary group U{2) which is not linearly 
related to the original one [21 . 

To be able to make sense of some of these nonlinear manipulations, it is conve- 
nient to consider a geometrical formulation of Quantum Mechanics. 

4. Geometric Quantum Mechanics 

The aim of this section is to present a brief summary of the set of the geometrical 
tools which characterize the description of Quantum Mechanics. For further details 
we refer the reader to El EH H El [HI E] . 

4.1. The states. To introduce the real manifold point of view, we start by replac- 
ing the Hilbert space Ti. with its realification TYr. In this realification process the 
complex structure on TL will be represented by a tensor J on TYr as we will see. 
The natural identification is then provided by 

Under this transformation, the Hermitian product becomes, for -0^,?/;^ G 7i 

To consider Tig just as a real differential manifold, the algebraic structures avail- 
able on ?i must be converted into tensor fields on Hr. To this end we introduce the 
tangent bundle TTYr and its dual the cotangent bundle T*7iR. The linear structure 
available in Hm is encoded in the vector field A 
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We can consider the Hermitian structure on TLn as an Hermitian tensor on TTi^. 
With every vector we can associate a vector field 

: Wm ^ TWr 4> {(j), i}) 

These vector fields are the infinitesimal generators of the vector group Hm. acting 
on itself. 

Therefore, the Hermitian tensor, denoted in the same way as the scalar product 
will be 

)((/)) = (V'l,V'2), 

i.e. the tensor field we are defining is translationally invariant. 

On the "real manifold" the Hermitian scalar product may be written as (V'l, ^2) = 
g{X^^,X^^) + icu{X^^,X^2), where g is now a symmetric tensor and u) a skew- 
symmetric one. It is also possible to write the associated quadratic form as a 
pull-back by means of the dilation vector field A extended as a map in an obvious 
way Hr ^ THu 8 TWr: 

The properties of the Hermitian product ensure that: 

• the symmetric tensor is positive definite and non-degenerate, and hence 
defines a Riemannian structure on the real vector manifold. 

• the skew-symmetric tensor is also non degenerate, and is closed with respect 
to the natural differential structure of the vector space. Hence, the tensor 

is a symplectic form. 

As the inner product is scsquilinear, it satisfies 

{ipi,iip2) = i{ipi,tp2), {■iipi,ip2) = -i{ipi,ip2)- 

This implies 

g{X^j^,Xji;^) = uj{JX^^,X^^). 

We also have that = —I, and hence that the triple {J,g,ui) defines a Kahler 
structure. This implies, among other things, that the tensor J generates both 
finite and infinitesimal transformations which axe orthogonal and symplectic. 

Linear transformations on the vector space Hr are converted into (1, l)-tensor 
fields by setting Ta where 

Ta : THm ^ THm cf>) ^ (V, Acj>). 

The association A — > T4 is an associative algebra isomorphism. The Lie algebra of 
vector fields obtains by setting Xa = Ta{A). Complex linear transformations will 
be represented by (1, l)-tensor fields commuting with J, or by vector fields Y such 
that LyJ = 0- 

For finite dimensional Hilbert spaces it may be convenient to introduce adapted 
coordinates on H and Hr. An orthonormal basis {le^)} of the Hilbert space allows 
us to identify this product with the canonical Hermitian product on C": 

(V'i,V'2) = ^('^i'i,efe)(efe, V2) 

k 

The group of unitary transformations on H becomes identified with the group 
U{n,C), its Lie algebra u{H) with u(n, C) and so on. 
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The choice of the basis also allows us to introduce coordinates for the realified 
structure: 

and write the geometrical structures introduced above as: 

J ^ dp^® dqk - dq^ (g) dpk g ^ dqk® dqk + dpk 'S) dpk w = dqk A dpk 

If we combine them into complex coordinates we can write the Hermitian struc- 
ture by means of z„ = qn + ipn- 

h = dzk (Xi dzk 
fc 

In an analogous way we can consider a contravariant version of these tensors. It 
is also possible to build it by using the isomorphism THw ^ T*Hr associated to 
the Riemannian tensor g. The result in both cases is a Kahler structure for the 
dual vector space Ti.^ with the dual complex structure J*, a Riemannian tensor G 
and a (symplectic) Poisson tensor fl: The coordinate expressions with respect to 
the natural basis are: 

• the Riemannian structure G — 'J2k=i ^ ^ ' ^ 



dqk dqk dpk dpk 

Sfc=l (, 

while the complex structure has the form 



the Poisson tensor fl — y^,',^ , ( A 



^ / A ® dqk - ^ ® dp,, 
\dpk dqk 

4.2. The observables. The space of observables (i.e. of self-adjoint operators 
acting on Ti) may be identified with the dual u*(7i) of the real Lie algebra u(7i), 
according to the pairing between the unitary Lie algebra and its dual given by 

A{T) = ^Tr^T 

Under the previous isomorphism, u*(7i) becomes a Lie algebra with product 
defined by 

[AB]^^{AB^BA) 
We can also consider the Jordan product: 

[A, B]+ = 2{A oB) = AB + BA 
with the associative product given by 

AB^ AoB+ -[A,BV 
n 

and defined on the complexified space 0[(7i). 

Both structures are compatible in the sense that the Jordan product is invariant 
under the linear transformations on u*(7i) defined by the adjoint action (of the 
Poisson product). As a result, u*{H) becomes a Lie- Jordan algebra (see [UlEl]). 

We can also define a suitable scalar product on Hermitian operators, given by: 

{A,B) = ^TrAB 

which turns the real vector space into a real Hilbert space. This scalar product is 
the restriction of the one on gl{n) defined as {M,N) = ^TrM^A^. 
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Moreover this scalar product is compatible with the Lie- Jordan structure in the 
following sense: 

([Ae],-B)u.(H) = (A[e,S])u.(H) + = (A[e,S]+)u.(W) 

These algebraic structures may be given a tensorial formulation in terms of the 
association A i-^ Ta on the tangent bundle Tu*{H). However we can also associate 
complex valued functions with linear operators A e by means of the scalar 

product 

gl{n)3A^U = l{iJ,Ai,)n. 
In more intrinsic terms we may write 

fA=l{9{A,XA) + MA,XA)). 

Hermitian operators give rise thus to quadratic real valued func;tions. 

The association of operators with quadratic functions allows also to recover the 
product structures on u{H) and u*{H) by means of appropriate (0, 2)-tensors on 
Hr. Indeed, by using the contravariant form of the Hermitian tensor G + given 
by: 

^.^,9 d d d d d . d d 
dzk ozk OQk oqk opk opk oqk opk 
it is possible to define a bracket 

{f.h}n = {f, h}g + t{f,h}^ 

We remark that and are not to be considered as derivations with respect 
to the complex coordinates introduced above but as complex valued smooth vector 
fields. 

In particular, for quadratic real valued functions we have 

{/a, .fB}g = /aB+BA — 2/aoB {/a, ~ —ifAB-BA 

The imaginary part, i.e. {•, - joj, defines a Poisson bracket on the space of func- 
tions. Both brackets define a tensorial version of the Lie-Jordan algebra for the set 
of operators. 

For Hermitian operators we recover previously defined vector fields; 

grad/^ = A; Ham/^ = iA 
where the vector fields associated with operators, we recall, are defined by: 

A-.Hr^THr 

iA-.Hf.^THR ^^iij,JA^P) 

We can also consider the algebraic structure associated to the full bracket {•,•}«, 
as we associated above the Jordan product and the commutator of operators to the 
brackets {•, ■}g and {•, respectively. It is simple to see that it corresponds to the 
associative product of the set of operators, i.e. 

{/a, /b}-h = {/a, ./elg + i{.fA, fBjuj = Iab+ba + Iab-ba = '^/ab 

This particular bilinear product on quadratic functions may be used to define a 
star product on quadratic functions 

Ua, fB}n = '^fAB = {dfA, dfB)H' = fA*fB 
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The set of quadratic functions endowed with such a product turns out to be a 
C*-algebra. 

Summarizing, we can reconstruct all the information of the algebra of operators 
starting only with real- valued functions defined on TYr. We have thus 

Proposition 1 ([H]). The Hamiltonian vector field Xf (defined as Xf — Cl{df)) 
is a Killing vector field for the Riemannian tensor G if and only if f is a quadratic 
function associated with an Hermitian operator A, i.e. there exists A = such 
that f ^ fA- 

Finally, we can consider the problem of how to recover the eigenvalues and 
eigenvectors of the operators at the level of the functions of Hr. We consider the 
expectation value functions associated to the operators as: 

Then, 

• eigenvectors correspond to the critical points of functions ca, i-C. 

deAi-ipa) = 

iff tpa is an eigenvector of A. We notice that the invariance of ba under 
multiplication by a phase U{1) implies that critical points form a circle on 
the sphere of normalized vectors if the eigenvalue is not degenerate. 

• the corresponding eigenvalue is recovered by the value eA(V'a) 

Thus we can conclude that the Kahler manifold {TL^, J,uj,g) contains all the 
information of the usual formulation of Quantum Mechanics on a complex Hilbert 
space. 

Up to now we have concentrated our attention on states and observables. If we 
consider observables as generators of transformations, i.e. we consider the Hamil- 
tonian flows associated to the corresponding functions, the invariance of the tensor 
G implies that the evolution is actually unitary. It is, therefore, natural, to consider 
the action of the unitary group on the realification of the complex vector space. 

4.3. The momentum map: geometrical structures on u*. The unitary action 
of U(TL) on TL induces a symplectic action on the symplectic manifold (TYr,^). By 
using the association 

F : He X u(W) ^ M {^,A) ^ \{iPMi>) = /»a(V'), 
we find, with Fa — fiA '■ '^r ^ I^, that 

{F{A),F{B)}^^^F{[AB]). 

Thus if we fix t/i, we have a mapping F{'4>) : u(7Y) — ^ M. With any element ijj ^ H 
we associate an element in u*(7i). The previous map defines a momentum map 

which provides us with a symplectic realization of the natural Poisson manifold 
structure available in u*(?i). We can write the momentum map from TYr to u*(7i) 
as 
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If we make the convention that the dual u*(7i) of the (real) Lie algebra u(7Y) 
is identified with Hermitian operators by means of a scalar product, the product 
pairing between Hermitian operators A £ u*{H) and the anti-Hermitian element 
T G u{H) will be given by 

A{T) = ^Tr(Ar) 

If we denote the linear function on u* (H) associated with the element iA e u(H) 
by yi, we have 

The pullback of linear functions on u* (Ti.) is given by the quadratic functions on 
Hm. associated with the corresponding Hermitian operators. 

It is possible to show that the contravariant tensor fields on Hm associated with 
the Hermitian structure are /^-related with a complex tensor on u*(H); 

/z*(G + iri) = R + iA; 

where the two new tensors R and A are defined by 

[R{A,Bm = (e,[AS]+>u* = iMaAB + BA)) 

and 

[A{A,B)m = [A,B]_U, = l-TrmB - BA)) 

Clearly, 

G{fi*A, ^i*B) + in{n*A, n*B) = ii*{R{A, B) + iA(i, B)). 
4.4. Example. Let us consider the manifold K^" and the dynamical system 

We can find alternative Hamiltonian structures connected to the factorization of 

0-i^M/7|-\ d o d 




-1 \ / 1 
\ I I 1 j i 



where 



1 ; - W -/.M / J, 

In the case of n = 2, by redefining position and momenta as {q,p) {q',p), 
according to the symmetric part of the decomposition we have: 



Let us now consider in the Hamiltonian function 

and the corresponding flow. 

Therefore, H = implies Pi + qi+P2 + Q2 = 1) i-^- a point in the unit sphere. 
A family of constants of the motion is provided by the set 

piq2-qiP2 PiP2 + qiq2 (pl + ql) - (pl + ql)- 
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It is simple to check that their Poisson bracket generate the Lie algebra of SU (2) 
and along with the dynamics they generate the Lie algebra of the full unitary group 
U{2). A foliation by two dimensional tori is given by the level sets of 

2,2 2,2 
Pl+Ql, P2 + l2^ 

on the open dense submanifold where they are independent. 

We notice that alternative Hamiltonian descriptions for F may be related to the 
immersion of F into different Lie groups or Lie algebras. For instance, we have seen 
that we may generate U{2) when we use the Poisson bracket defined as 

{qa,Pb} = Sab {qa,qb} ^ ^ {Pa,Pb}- 

But if instead we consider 

= {'Z2,j32} = -1, {qa,qb} {pa,Pb}; 

the same constants of the motion will generate the Lie algebra of the group SU{1, 1). 
This would correspond to the case of the factorization of the matrix A associated 
with the linear dynamics into a skew-symmetric and a symmetric one which is not 
degenerate but it is not positive definite. 

It is also possible to consider the dynamics as part of an infinite dimensional Lie 
algebra (deformations) . 

In general we look for factorizations of the vector field by means of a Poisson 
tensor and a differential one form (F — A{dH)), in this way both A and dH need 
not be associated with matrices. In particular if we consider the (1, l)-tensor field 

with any constant of the motion of the form 

F = F{pI + qf,pl + ql,piq2 - qiP2,PiP2 + 9i'?2), 
we can associate an invariant one form |47j 

dp = dj'F, 

By a proper choice of the function F, we can make the two form ddxF non degen- 
erate, i.e. we can choose F in such a way that the two form 

LUp = ddxF, 

is an invariant symplectic two form. 
Moreover 

LrdxF = irddpF + dirdxF = 0, 

implies that 

irddrF = —d{irdTF), 

with 

i^drF = dF{T{T)) = -iydF{A) = -i^LaF 

because T(F) = —vA, A being the dilation vector field. 

Any other invariant (1, l)-tensor field will give rise to alternative symplectic 
structure with the same procedure. 

More generally, if Xj and A"j+„ are pairwise commuting vector fields [Xj, Xk] — 
Vj, k — I, - ■ ■ , 2n, with 

LrXj — Xjj^n L-pXj^n — —Xj] 
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the one form 

3 

is F-invariant and under proper choices of F provides a syniplectic structure 

(jjp = dOp. 

In addition, if we set 

we find a new Hnear structure and a new Kahler structure invariant under T: 

LvQj = Pj LrPj = -Qj. 
The new Haniihonian function is given by 

\{{Lx,FY + {Lx,,,Ff) = \Y.iQ'3+Pl)^ 

3 

for any constant of the motion F ^ which makes LOp not degenerate. 

This construction of ahernative Kahlerian structures applies therefore to the 
case of quantum systems along the lines of [21] • 

By using the Lie algebra generated by 

we can consider a larger algebra by using an infinitesimal generator of an element 
in the automorphism group of the algebra itself 

3 

In general, any quadratic Hamiltonian may be considered as an element of the 
automorphism group, the inhomogeneous symplectic group. In particular we may 
find subalgebras or subgroups which contain just one quadratic Hamiltonian and 
the full linear part. This would be a generalization of what is known as the oscillator 
group, which is four dimensional |52j . 

It is interesting to notice that the oscillator group allows for a "contraction" by 
letting J/ ^ and therefore it can reproduce the free particle case. On the other 
hand, considering the unfolding procedure described in HJ |3] it is also possible to 
describe the Hydrogen atom (Kepler problem) by using a nonlinear Harmonic oscil- 
lator (i.e. by introducing an energy-dependent period to be able to accommodate 
for the energy-period theorem when dealing with the central- force problem). 

5. The notion of independence 

5.1. Proposal based on the geometric framework. Let us go back again to 
the problem discussed above originated by the von Neumann theorem, but from 
a more constructive perspective. Let us see how the tools we introduced within 
the geometric description of quantum mechanics allow us to introduce a notion of 
independence which may be convenient for our purposes. 

Let us consider an even dimension vector space V and a linear differential equa- 
tion associated to an operator A: 

i, — Ax X 

A is defining a vector field on V . We can consider: 
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• there is a symplectic form ujonV such that the vector field is Hamiltonian. 

• Solve the inverse problem of Hamiltonian mechanics to find invariant sym- 
plectic forms 

If V is finite dimensional, we can obtain metric tensors (by using the Hessian 
of the Hamiltonian function) such that the vector field is a Killing vector field for 
g. By proper rescaling the metric tensor defined in this way, we may combine co 
and the rescaled g to define an Hermitian product on V which is preserved by the 
Hamiltonian vector field. This allows us to introduce a set of quadratic functions 
on V: 

B ^ fB = {x,Bx), 

associated with Hermitian matrices with respect to the newly defined Hermitian 
products. 

The set of (quadratic) functions can be endowed with different algebraic struc- 
tures: 

• the pointwise product 

• the Poisson bracket associated to the symplectic form 

• the non-local product which translates the product of operators at the level 
of quadratic functions: 

AB^fA^fB 

• the skew-symmetric part of the above, which encodes the commutator al- 
gebra of the operators: 

Let us now discuss how to introduce a concept of functional independence. We 
have two possibilities: 

• to consider the independence in the sense of the pointwise product and to 
use the standard differential calculus. Two functions are then said to be 
independent iff 

dfA AdfB^O 

• to consider independence with respect to the ★ product: in that case a non- 
commutative differential calculus is required and the issue becomes much 
more involved (see [5TJ [SUl IHl Si] ) . 

In the finite dimensional case it is simple to see that the notion associated with 
the pointwise product is convenient and provides a reasonable framework. Besides 
it has the advantage of being easily extensible to the case of infinite dimensional 
systems, once: 

• a metric tensor compatible with the symplectic form has been defined. The 
definition can be taken from |40j . where the operator A is used: once the 
symplectic form has been chosen the Hamiltonian quadratic function /a is 
considered. This is a quadratic function, and hence it can be used to define 
a metric tensor 

Then, there are two possibilities: 

— it is positive definite: in this case it is used directly to write the metric 
tensor as 

g = [ dq dp ^ Hess f a 
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— if it is not positive definite it can be proved (|42]) that the space V can 
be decomposed into subspaces (say Vi and V2) so that the functions 
may be written as 

Ja = Iai — !a-2 

where A\ is a matrix with entries only on Y\ and A2 with entries only 
on V2. A pseudo- metric (non definite metric) is thus obtained. The 
corresponding Hermitian structure will be pseudo-Hermitian on the 
total Hilbert space, but integrability may be then discussed on each 
Hilbert subspace where the metric may be considered to be Hermitian. 
• the metric tensor allows to define suitable functions corresponding to the 
powers of the operator A. This is not a trivial issue, since depending on 
the domain of A, the powers need not have an associated function (the 
corresponding quadratic expression need not be convergent). Hence the 
freedom in the choice of the metric may be used to force all powers of the 
operator A to be trace-class. 

5.2. An example. Let us consider an application of this construction to a partic- 
ular example. Consider a finite level system defined on H and an operator A and 
its square . It seems obvious to consider both operators as dependent, from the 
point of view of the associative operator product. But let us consider their inde- 
pendence from the point of view of the associated quadratic functions with respect 
to the differential calculus associated with their commutative product. Consider 
therefore 

jA = (V-I^l^) !a- = (^l^'l^) 

If the operator A is self-adjoint with simple eigenvalues, we can consider the 
basis of eigenstates, and write: 

/a = XI ^^"^ = 5Z ^1 

3 j 

If we compute the differential of each function: 

dfA = ^j{zjdzj + Zjdzj) dfA^ = ^'j{zjdzj + Zjdzj) 
j 3 
And these two one forms satisfy 

dfA A dfA^ = X^^^J'^fc ^ \]\k)zjZkdzj A dzk) = ^ \j\k{\k - \j)zjZkdzj A dz^) ^ 

jk jk 

Thus we see that, in the generic case, these two functions can be considered to 
be functionally independent. The condition for this, in general, will depend on the 
spectrum. 

Remark 

We try to explain in a nutshell the main ideas involved in previous constructions. 

Any algebra A may be considered as a vector space V along with a binary 
product B : V X V ^ V. Given any linear vector space V, we may "geometrize" it 
by considering the embedding V ^ Lin{V* ,M.) C T{V*), i.e. any element u of V 
is thought of as a linear function on the dual space V* . If a G , we have 

u{a) — a(w), 
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i.e. u is mapped into u. Once we identify u with a function on V* , it makes sense 
to consider polynomial functions and construct more general functions by means 
of the Weierstrass theorem. On this space of functions we may develop a standard 
exterior differential calculus. 

On the other hand, the binary bilinear product available on V , allows to define 
a product on V, by setting 

{v-ku){a) = a{B{v,u)). 

This product is non-local and would require a non-commutative exterior differential 
calculus. 

To clarify our statements, let us consider the algebra of n x n matrices M(n,M). 
We may consider its dual space generated by 

{ejk = \ej){ek\}, 

constructed out of an orthonormal basis in some n-dimensional vector space E. 
Clearly, 

ejkiA) = {ek,Aej) 
defines the dual pairing. We may define thus 

and the products 

{A-B){ejk) = ejk{A)-ejk{B) 

and 

{A*B){ejk) = ejkiA-B). 

The first one is the usual Hadamard (or Shur) product among matrices and is 
commutative and local. The exterior differential calculus associated to it is the 

2 

standard one on M" . 

The second product is the standard row-by-column product, and is non-local 
and non-commutative. 

Thus, in a sense, the first product is originated only from the vector space 
structure of the algebra, while the second one is specific of the binary bilinear 
product which defines it. Our suggestion is thus to use the exterior differential 
calculus associated with the first (local and commutative) product to define the 
notion of functional independence of observables in Quantum Mechanics. 

6. Examples: Harmonic oscillators and their deformations 

From our previous discussion on classical integrability, it is obvious that Har- 
monic oscillators of any dimension will provide examples of integrable systems, also 
in the quantum framework. We will use them as a reference model to which all other 
quantum integrable systems should be related (those having a discrete spectrum, 
of course). 

Now we are going to review some properties of the system and we will present 

then two examples of quantum integrable systems which exhibit the properties we 
have discussed above: the deformed Harmonic oscillator and the Coulomb problem 
in two and three dimensions. 
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6.1. The standard Harmonic oscillator. As a preliminary simple example, we 
consider the Harmonic oscillator. The Hamiltonian is well known: 

supposing h — \ — uj. 

The dynamics is linear and can be integrated easily. This system is obviously 
integrable from our point of view, and in both schemes (vector fields or operator 
subalgebras) the mapping <i> is the identity. 

To deal with the Harmonic oscillator one introduces creation and annihilation 
operators a, and the vacuum state |0) and builds the states of the Fock space as 



'71! 

These are eigenstates of the number operator n = a) a, and thus define the eigen- 
states of the Hamiltonian which is written as 

The operators satisfy the commutation relation 

[a, at] =1 

and the number operator (n = a^a) thus generating together with the identity, the 
well known Heisenberg-Weyl group {H4). 

The Hilbert space structure is fixed by the usual scalar product 

We notice that in the geometric framework these two objects determine the 
Riemannian and Poisson structures of our description. 

In the Heisenberg picture the equations of the motion are written as 

a — —ia a) = ia) 

It is well known that the number operator classifies the irreducible representa- 
tions, as the Casimir corresponding to the subgroup C/(l) C -ff4. The corresponding 
Hilbert space is the well known Fock space V . 

Our aim is to show how the same dynamics can be described on the space of 
states endowed with different Hermitian structures (therefore on different Hilbert 
spaces). The two different structures are related by a nonlinear transformation 
which shows how these types of transformations arise within the quantum frame- 
work as described in |37l |MI . 



6.2. A non trivial combination: the /— deformed Harmonic oscillator. Let 

us consider a construction similar to the previous case, but with some novelties. 
Let us introduce the operators 

A = aJ{h) A^ = J{n)a\ 

where / is a certain smooth function of the number operator (or equivalently, a 
function of the energy of the system) . As the number operator is a constant of the 
motion for the dynamics of the Harmonic oscillator, it is clear that the dynamics 
in this new operator basis has the same form: 

A ^ -lA it ^ 
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Associated with these operators we can also define a set of vector states (the 
vacuum state |0) is the same as before): 

\N)=^-yi^\0), 



which are the eigenstates of a number operator defined as 

N = A^A. 

A scalar product on this new set is defined also as 

{N\M) = 5nm. 

Wc sec thus that there is a one-to-one correspondence between the two sets of 
states, since the mapping 

T\:\k)^\K) k,Ke No 

is clearly bijective. It is important to remark, though, that the linear structures are 
not preserved by this mapping, which is nonlinear, i.e. n(|ni) + 1^2)) 7^ n(|ni)) + 
n(|n2)). We shall come back to this point in the next section. 

We have thus two difi'erent Hilbert space structures associated with the initial 
system. We have here a situation very similar to the classical one: the "abstract" 
Heisenberg-Wcyl algebra is realized in difl^crcnt manners and gives rise to algebras of 
operators related by an intertwining operator which is not a linear transformation. 

The situation can be summarized as follows: 

• With respect to the initial Hilbert space structure, we can write the com- 
mutation relation of a and as 

(n, [a,a^]m) = Snm- 

An analogous relation exists for the operators A and A^ with respect to the 
new Hermitian structure: 

{N, [A,A^]M) = 

• But if we consider A and as operators acting on the original Fock space, 
the commutation relations with respect to this action have the following 
form 

(n, [A, At]m) = ((n + l)fin + 1) - nfinf)5nm. 
We can write thus that, with respect to the original Fock space structure, 
[A,A'<] = {n+l)f{n+l)-nf{nfl 

However the equations of motion of the Harmonic oscillator, i.e. the vector field, 
are the same in both cases, but the choice of one Hilbert space structure or the 
other makes the description quite different when written in the "wrong" framework. 
And notice also that both Hilbert space structures arc related to each other by a 
transformation which is nonlinear. The corresponding transformations at the level 
of the symplectic structures (i.e. the transformations which relate the commutation 
relations) define a non-canonical nonlinear transformation. 

We can give some more details of the construction if we choose a particular 
example of the mapping /: 
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6.2.1. One dimensional example. A typical example of this construction is provided 
by the so-called g-deformed Harmonic oscillators, defined by the expressions: 

sinh nh ^ 

sinh h 

where n is the undeformed number operator introduced above. Analogously the 
creation-annihilation operators are written in terms of the undeformed ones: 

aq = af{n) a+ = /(n)a+ 

where 

Notice that this transformation is invertible, since we can write 



a = A, 



log (iV sinh sinh h + 1 



But it is simple to read from here the non-linearity of the mapping 11 ; |fc) ^ \K) 
introduced above, since the relation for the number operator is thus written as: 



1 



n = ^ log sinh H + ^T^sinh^ H+l 

The mapping 11 is just read from here since it is defined on the set of eigenstates 
of these two operators (n and TV'). 

The expression for the new commutation relation in terms of the old Hilbert 
space structure corresponds to 



,2 



[A,A'<] = (^N{coshn- 1) + \In^ sinh^^i-hl^ I 
For this scheme, we define the Hamiltonian of the g-deformed oscillator to be: 

Hq = ^log ^7Vsinhn-h xj N"^ sinh^ft+l-h . 

With this operator and the original Hilbert space structure, we can write the asso- 
ciated vector field: 

[A,H,] = A, 

which gives the original dynamics written in terms of the new operators. 

6.2.2. Two dimensional case. For the cases with more than one oscillator we have 
the option of keeping or not the U{N) symmetry of the undeformed problem. If 
we prefer to keep the symmetry, we define the deformed creation and annihilation 
operators as: 

aq = af{n) a+ = f{n)a~^ 

bq = a fin) h+ = f{n)a+ 

where now n = Ua + ni, (to break it, one should choose for instance only one of the 
number operators for each deformation). 
The Hamiltonian is now obtained as: 

H=hnf{n) + {n+2)f{n + l)] 
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which exhibits the kind of combination of models we have exposed above. It is 
a combination of Harmonic oscillators with a function which is a "constant of the 
motion" on each description. 

If we want to break the symmetry and define a different combination for the 
Hamiltonian, we define the creation and annihilation operators as: 

aq = af{na) a+ = f{na)a^ 
bq = afim) 6+ = I{nb)a^ 
in such a way that the Hamiltonian becomes: 

H = ^K/2(n„) + nj,f{jn,) + (n, + l)/^^ + 1) + (^fc + 1)/'K + 1) 

In this case is even more clear the way models are combined. We have two 
different oscillators, and the construction of a Hamiltonian with different factor for 
each representation of them. 

The mapping <i> relating the two descriptions, acts only on each particular repre- 
sentation of the algebra of each particular oscillator. It transforms the algebra 
and leaves invariant the Hilbert space, which is constructed as the product of both 
Fock spaces in the usual way. The global transformation is constructed as: 

where is the one dimensional transformation and it is equal to the identity. 

6.3. The Coulomb problem. As an example of the potential applications of the 
composition of models to the study of quantum problems let us consider a simple 
example in the Schrodinger picture. 

It is well known (see for instance [5U]) that the system may be described with 
the help of a combination of Harmonic oscillators for the bound states in the D = 2 
and D = 3 cases. From our point of view this implies that the system is integrable 
at the quantum level (for these dimensions and regimes). 

For both cases we decompose the Hilbert space in the eigenspaces of the Hamil- 
tonian, let us call them Ue^ ■ Now, on each subspace we perform a change of 
coordinates that maps Ue^ on a suitable subspace of a Fock space.: 

• For the two dimensional case the change of coordinates is simply the trans- 
formation into parabolic cylinder coordinates: 

X — 2^"^^ ~ y — uv 

where the point (— w, —v) and (u, v) are identified. This change transform 
the eigenvalue equation 

1 / d'^ip d'^ip\ k 



into 



^2j, Pt2, 



■0 = Ei\) 



1 1 2. 2,, p,. 



2 

where the wave-functions should verify '0(— u, — w) = ?/'(u,v), where uP' — 
2\E\ = || and E' — k = ruv where n has to be odd to satisfy the parity 
condition of the wave functions. We construct thus a mapping from the 
Hilbert space of bound states of the two dimensional Coulomb problem 
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into the subspace of the product Fock space Vu (E> Vv corresponding to an 
odd number of particles. There is a dependence of the mapping on the 
Hamiltonian, which produces a different frequency for each oscillator. 

■■ Ue^ IC''' c K K (19) 
The complete mapping is the union of these e which are formally identical: 

• For the three dimensional case, the situation is analogous (see |20j as the 
original reference) . Now we make use of the su(2)©su(2) symmetry algebra 
of the Coulomb problem (corresponding to suitable combinations of the 
angular momentum and the Runge-Lenz vector if we define it in terms 
of generators) . We construct then a two dimensional Harmonic oscillator 
associated with each of the su(2) group, thus obtaining a four dimensional 
isotropic Harmonic oscillator. 

The change of coordinates is constructed in two steps: 

Xi = fipcoscj) X2 = fivsincj) = 2^1-''^ ~ ^'^) 
= ^ cos rj^j_ = ^. sin 4> 
— V cos 4> Vv = V sin 

where again = 2\E\ (the change is different for each eigenspace of the 
Hamiltonian). These changes transform now the system (this time consid- 
ered in three dimensions) into: 

where 

and E' = k = Lu{nf^ + 11^ + 1). 

The mapping should be understood now applying each eigenspace of the 
Hamiltonian of the hydrogen atom, representation by representation of the 
angular momentum and the Runge-Lenz vector, onto the corresponding 
representations of the Hamiltonian of the Harmonic oscillator. 

We can write analogously to the previous case: 

: Ue^ ^ V;°f C ® K (20) 

where now Va is the Fock space corresponding to a two dimensional oscil- 
lator. For the whole Hilbert space we can write again: 

A different approach, where starting from the classical situation the construction 
of the quantum one is directly related to the geometry of the Lagrangian description 
is provided in [S]. 
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7. Conclusions and outlook 

In this paper we have pointed out various problems that a definition of quantum 
integrabihty has to face. The main one seems to be the definition of functional 
independence which, in our present understanding, requires a differential calculus 
on commutative algebras. This calculus deals with derivations and is fully captured 
by the notion of differential manifold, be it finite or infinite dimensional. For this 
reason we believe that the geometrization of quantum mechanics can be a good 
starting point. In this approach we find that alternative products on functions 
coexist: 

• a non-local product (on quadratic functions) which captures the essence 
of quantum mechanics in terms of indetermination relations and a rule for 
composing systems, 

• a local commutative product which allows for a differential calculus in terms 
of derivations which are local. 

At the moment, it seems that the latter is useful, if not unavoidable, to state 
functional independence of "observations" (quadratic function version of the ob- 
servables) . 

Usual treatments of integrability ( |49] ) at the quantum level deal with Hilbert 
spaces considered as spaces of square integrable functions on some "configuration 
space" . It seems that the notions we are advocating are provided by Lie groups and 
their unitary representations, along with the associated representations of their Lie 
algebras (and their enveloping algebras) on the set of analytic or smooth vectors 

[mile]. 

The convolution product and the pointwise product available on the group seem 
to be related to the two products we are identifying. All these considerations 
appear very neatly in the so called Quantum Mechanics on phase space (the Wigner- 
Weyl formalism, see [S]) There the phase space is an Abelian vector group and the 
convolution product (Moyal product) is associated with its central extension defined 
with the help of a symplectic structure). The algebra of differential operators is 
identified as the algebra of operators acting on square integrable functions defined 
on Lagrangian subspaces of this symplectic vector space. 

The representation of the enveloping algebra allows us to deal with suitable 
deformations (quantum groups and Yang Baxter relations). 

The transformation from one Lagrangian subspace to another one requires the 
introduction of pseudo-differential operators. Nonlinearly related actions of E^" 
on the phase space, when considered as Weyl systems [2T], give rise to nonlinear 
transformations in the quantum framework. 

Some of these statements need further elaboration to be able to convert them 
into propositions and theorems. We hope to come back to these questions in the 
near future. 

Another spin-off of this geometrization procedure is the possibility of using meth- 
ods and conceptual constructions elaborated in the framework of classical physics 
(vector fields, Poisson brackets, Riemannian vector fields) also in the framework of 
quantum physics, with the additional intervention of additional structures which 
may be traced back to the existence of the fundamental constant of nature provided 
by Planck's constant 
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The difficulties we have to face can be traced back to the lack of a mature 
mathematical treatment of infinite-dimensional differential manifold even though 
very good books are already available ( |33i ITUl [50] ) . Our hope is to be able to 
apply in the very near future these methods to study quantum control problems 
using as much as possible the experience available in the "classical framework" . 
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